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ABSTRACT 


Detailed  measurements  of  the  turbulent  multiphase  flow  associated  with  wave 
breaking  present  a  unique  instrumentation  challenge.  Measurement  systems  must  be 
capable  of  high  sampling  rates,  large  dynamic  ranges,  as  well  as  be  capable  of  making 
measurements  in  water,  air  and  optically  opaque  regions.  An  experiment  was  performed 
on  Carriage  2  in  the  Deep  Water  Basin  at  the  Naval  Surface  Warfare  Center,  Carderock 
Division,  (NSWCCD)  in  October  2007  to  measure  various  characteristics  of  the  breaking 
wave  generated  from  a  submerged  ship  transom.  The  primary  objective  of  this  work  was 
to  obtain  full-scale  qualitative  and  quantitative  flow  field  data  of  a  large  breaking  transom 
wave  over  a  range  of  conditions,  specifically  transom  drafts  and  Froude  numbers. 

Several  types  of  measurements  were  made  of  the  transom  stem  wave.  Sinkage 
and  trim  were  measured  using  two  string  potentiometers.  Drag,  vertical  and  side  forces 
were  measured  using  block  gages.  To  quantify  the  spray  and  free  surface  deformation, 
several  techniques  were  used,  including  a  scanning  LiDAR  system,  laser  sheet  flow 
visualization  (Quantitative  visualization  or  QViz),  and  Senix  Ultrasonic  acoustic  distance 
sensors.  Additional  measurements  were  made  using  the  Nortek  Acoustic  Wave  and 
Current  (AWAC)  profiler,  which  measured  velocity  and  acoustic  backscattcr  through  the 
water  column.  An  array  of  impedance  void  fraction  probes  was  also  used  to  measure  the 
entrained  air  at  various  locations  and  depths  behind  the  stem. 
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INTRODUCTION 


The  physics  of  the  transom  stem  wave  continues  to  be  of  great  importance  to 
understanding  both  ship  breaking  waves  and  bubble  wakes.  The  full-seale  breaking 
transom  stem  wave  is  a  complex  non-linear  turbulent  flow  field,  and  while 
Computational  Fluid  Dynamics  (CFD)  codes  have  demonstrated  improved  capability  in 
predicting  the  large-scale  Kelvin  wave  structure  for  a  variety  of  naval  craft,  the  ability  of 
CFD  codes  to  predict  the  short-scalc  surface  evolution  and  the  energy  dissipation 
involved  in  breaking  regions,  spray  sheets,  and  turbulence  has  not  yet  been  validated  and 
remains  a  challenge.  The  primary  objective  of  this  work  is  to  obtain  full-scale  qualitative 
and  quantitative  flow  field  data  of  a  large  breaking  transom  wave  over  a  range  of  Froude 
numbers. 

An  experiment  was  performed  in  the  Deep  Water  Basin  at  the  Naval  Surface 
Warfare  Center,  Cardcroek  Division,  (NSWCCD)  in  October  2007  to  measure  various 
characteristics  of  the  breaking  wave  generated  from  a  submerged  ship  transom.  Several 
types  of  measurements  were  made  of  the  transom  stem  wave.  Sinkage  and  trim  were  was 
measured  using  two  string  potentiometers.  Drag,  vertical  and  side  forces  were  measured 
using  block  gages.  To  quantify  the  spray  and  free  surface  deformation  several  techniques 
were  utilized,  including  a  scanning  LiDAR  system,  Digital  Light  Projection  (DLP) 
enhanced  Quantitative  Visualization  (QViz),  and  Senix  Ultrasonic  acoustic  distance 
sensors.  Additional  measurements  were  made  using  the  Nortck  Acoustic  Wave  and 
Current  (AWAC)  profiler,  which  measured  velocity  and  acoustic  backscatter  through  the 
water  eolumn.  An  array  of  impedanee  void  fraction  probes  was  also  used  to  measure  the 
entrained  air  at  various  locations  behind  the  stem  and  various  water  depths. 


EXPERIMENTAL  APPROACH 


Model  Description  and  Facilities 

This  experiment  was  performed  by  towing  Model  5673  in  the  Deep  Water  Towing 
Basin  on  Carriage  2.  The  basin  is  approximately  22  feet  deep,  1886  feet  long  and  50.96 
feet  wide,  with  a  maximum  carriage  speed  of  33.8  ft/s  (20  knots)  (Saunders,  1).  The 
model  was  towed  using  a  tow  post  located  270  in  (6.9  m)  forward  of  the  aft 
perpendicular.  A  grasshopper  was  used  90  in  (2.3  m)  forward  of  the  aft  perpendicular  to 
fix  the  model  in  yaw,  while  still  allowing  it  to  pitch. 

Model  5673  (shown  in  Figure  1,  Figure  2  and  Figure  3)  has  a  transom  stem  and  was 
designed  to  minimize  the  generated  bow  wave  so  that  the  transom  wake  eould  be  more 
effectively  investigated.  The  model  is  about  30  feet  (9.1  m)  long,  with  a  maximum  beam 
of  5  feet  (1.5  m).  Figure  4  shows  a  plan  view  of  the  model  mounted  on  Carriage  2. 
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Figure  1.  Image  of  the  transom  model  geometry. 


Figure  2.  Plan  and  profile  views  of  the  transom  model  geometry. 
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NSWC-  Carriage  2 

30-ft  Model,  facing  east 


Figure  4.  Plan  view  of  model  mounted  under  Carriage  2. 


Test  Conditions 


Though  it  was  intended  to  test  the  model  at  a  fixed  draft  and  trim  position,  the 
large  forces  generated  prevented  this  towing  configuration.  Instead,  the  model  was  tested 
as  fixed  in  heave  and  free  to  pitch.  The  pitch  magnitude  was  recorded  during  testing. 
Four  different  speeds  were  tested,  including  5,  7,  8  and  9  knots.  Table  1  shows  the  test 
conditions  for  this  experiment,  along  with  the  length  and  draft  Froude  numbers,  where 
length  Froude  number  is  defined  as: 


v 


(1) 


and  draft  Froude  number  is  defined  as: 


FnD  = 


v 

7^ 


where  v=  model  velocity 
g=gravitational  acceleration 
L=length  of  model  (30  ft  for  this  model) 
D=draft  at  the  transom 


(2) 
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For  the  5  and  7  knot  conditions  the  transom  was  partially  wet,  and  for  the  8  and 
9  knot  conditions  the  transom  was  entirely  dry.  Literature  suggests  that  a  transom  stem 
vessel  will  experience  a  dry  stem  (fully  ventilated)  at  draft  Froude  numbers  above  2.5 
(Maki  et.  al.,  2,  and  Faltinsen,  3).  In  this  experiment,  the  transom  stem  is  dry  at  a  slightly 
lower  draft  Froude  number  of  about  2.1 1 . 

The  model  was  tested  at  five  different  longitudinal  positions  relative  to  the 
carriage  to  accommodate  the  various  data  collection  systems  that  were  used.  These 
positions  arc  listed  in  Table  2. 


Table  1.  Test  Conditions. 


Speed 

(knots) 

Length 

Froude  Number 
(FnL) 

Draft 

Froude  Number 
(Fn0) 

Transom 

Condition 

5 

0.27 

1.4 

wet 

7 

0.38 

1.9 

wet 

8 

0.43 

2.1 

dry 

9 

0.49 

2.3 

dry 

Table  2.  Model  positions  relative  to  carnage 


Position 

Distance  from  the  east  end  of  the  floating  girder  to  the 
east  end  of  the  tow  post  (along  the  girder,  inches) 

i 

44.25 

2 

63.375 

3 

9.125 

4 

21.625 

5 

30.375 

Instrumentation 


Standard  Video  and  Still  Imaging 

Four  standard  frame  rate  (30  fps)  video  cameras  were  used  to  record  video  during 
the  test.  One  camera  captured  a  bow  view,  and  the  other  three  captured  images  of  the 
free  surface  aft  of  the  model.  Also,  a  digital  still  camera  was  used  to  record  the  visual 
appearance  of  the  free  surface  aft  and  around  the  transom  model. 


Block  Gages 

Three  calibrated  4  in  (10  cm)  block  gages  were  used  under  the  tow  post  to 
measure  the  lift,  drag  and  side  forces.  An  additional  4  in  (10  cm)  block  gage  was  used 
under  the  grasshopper  to  monitor  side  forces.  The  block  gages  were  calibrated  by 
NSWCCD,  Code  5800,  following  standard  procedures. 
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String  Potentiometers 


Model  trim  was  measured  using  string  potentiometers  loeated  at  the  bow  and  stem 
of  the  model.  The  distance  between  the  string  potentiometers  was  216.8125  in  (5.5  m). 
The  forward  string  pot  was  loeated  288.125  in  (7.3  m)  forward  of  the  aft  perpendicular 
and  the  aft  string  pot  was  located  71.3125  in  (1.8  m)  forward  of  the  aft  perpendicular. 
The  string  potentiometers  were  calibrated  by  NSWCCD,  Code  5800,  following  standard 
procedures. 


LiDAR 


Light  Detection  and  Ranging,  or  LiDAR,  is  a  remote  sensing  system  used  to 
collect  topographic  data.  The  LiDAR  system  used  during  the  2007  Transom  Test 
contains  a  single  Riegl  pulsed  laser  scanner  (LMS-Q140-80i)  and  a  four-sided  mirror 
which  spun  to  deflect  the  laser  at  different  angles  along  a  single  line.  The  time  for  the 
reflected  pulse  to  echo  back  to  the  sensor  receiver  is  used  to  calculate  distance.  The 
range  accuracy  of  the  LMS-Q140-80i  unit  is  generally  +/-  0.8  in  (2  cm),  which  typically 
scans  in  a  +/-  40  degree  sweep  at  a  laser  pulse  frequency  of  30  kHz. 

LiDAR  data  was  collected  at  a  rate  of  20  Hz.  The  system  was  mounted  above  the 
transom  of  the  model  on  a  traverse  attached  to  the  carriage,  in  an  effort  to  measure  the 
surface  wave  field  generated  in  this  region.  The  LiDAR  system  was  mounted  12.7  ft 
(3.87  m)  above  the  still  waterline  of  the  tank  and  0.82  ft  (0.25  m)  starboard  of  Model 
5673.  Table  3  shows  the  locations  of  the  LiDAR  for  the  various  data  sets  that  were 
eollceted,  where  the  locations  are  relative  to  the  stem  of  the  model. 


Table  3.  LiDAR  locations  and  conditions  for  data  collection  (locations  are  relative  to  the  stem  of  the 

model). 


Aft  Location  (inches) 

Speed 

<w») 

59.8" 

67.8” 

68.8" 

76.8" 

77.3" 

82.3" 

85.8" 

86.3" 

90.3" 

943 

99.8" 

103.8" 

108.3" 

116.3" 

5 

✓ 

✓ 

✓ 

V 

7 

✓ 

✓ 

8 

✓ 

✓ 

✓ 

✓ 

✓ 

S 

✓ 

✓ 

9 

✓ 

✓ 

✓ 

S 

V 

h” — 

=  Dry  Transom 


DLP-Enhaneed  OViz 

A  non-intrusivc  optical  technique,  Digital  Light  Projection  (DLP)  enhanced 
Quantitative  Visualization  (QViz),  was  developed  to  pursue  free  surface  measurements  at 
high  spatial  and  temporal  resolution.  A  DLP  was  used  to  project  a  laser  light  sheet 
perpendicular  to  the  free  surface,  and  video  cameras  were  used  to  collect  digital  images 
of  the  intersection,  representing  instantaneous  eross-seetions  of  the  wave  shape.  The 


6 


latter  aspect  of  the  system  operation  was  similar  to  previous  versions  of  the  QViz  system, 
details  of  which  are  given  in  Furey  and  Fu  (4),  and  Rice  et  al.  (5).  However,  the 
introduction  of  DLP  technology  permitted  several  vertical  light  sheets  to  be  scanned 
throughout  a  test  run,  allowing  the  free  surface  to  be  effectively  mapped  over  a  desired 
area. 

The  novel  projection  optics  of  the  DLP-enhaneed  QViz  system  used  a  Digital 
Mieromirror  Device  (DMD),  an  optical  semiconductor  instrument.  The  DMD  device 
(Texas  Instruments,  DMD  Discovery  1100)  contains  an  array  of  1024  by  768 
mieromirrors.  In  the  system  configuration  used,  it  was  controlled  using  a  USB  interface 
to  project  lines  onto  the  free  surface  that  were  1024  mirror  pixels  long  and  approximately 
10  mirror  pixels  wide.  A  timing  signal  was  sent  from  the  DMD  device  to  the  video 
cameras  so  that  the  projected  images  were  synchronized  with  standard,  30  fps,  video 
cameras.  The  projection  optics  were  mounted  on  a  2.5  ft  by  2.5  ft  (0.75  by  0.75  m)  square 
optical  breadboard  which  was  required  to  be  located  directly  above  the  desired  wave 
region. 

Large  amounts  of  free  surface  image  data  were  collected  at  five  locations.  For 
each  location,  lines  were  projected  to  scan  a  measurement  area  of  approximately  1.3  by 
1 .6  ft  (0.4  by  0.5  m)  within  the  duration  of  each  run.  Approximately  60  eross-seetions 
were  obtained  for  each  line.  Three  locations  were  centered  transversely  at  0.66  ft  (0.2  m) 
starboard  of  the  model  centerline,  and  longitudinally  at  6.2,  6.9,  and  7.9  ft  (1.9,  2.1  and 
2.4  m)  aft  of  the  model  transom.  Two  additional  locations  were  centered  longitudinally 
at  5.9  ft  (1.8  m)  starboard  of  the  model  centerline,  6.6  and  8.2  ft  (2.0  and  2.5  m)  aft  of  the 
model.  These  two  measurement  regions  represented  an  effort  to  collect  data  across  the 
edge  of  the  breaking  region,  or  shoulder,  of  the  wake. 

Senix  Ultrasonic  Sensors 


Seven  Senix  Ultrasonic  sensors,  whieh  are  non-eontaet,  acoustic  instruments  for 
measuring  distances  through  air,  were  used  to  collect  longitudinal  waveeut  data.  A  truss 
section  (wave  boom)  cantilevered  from  the  basin  wall  over  the  water,  provided  a  structure 
on  whieh  the  sensors  were  mounted,  as  shown  in  Figure  5.  The  wave  boom  extends  22.4 
ft  (6.83  m)  from  the  basin  wall,  whieh  places  the  end  of  the  wave  boom  approximately  3 
ft  (0.91  m)  short  of  the  basin  centerline.  A  photosensor  was  set  to  trigger  data  collection 
when  the  forward  perpendicular  of  the  model  24.979  ft  (7.61  m)  from  the  sensors.  Wave 
elevation  data  was  collected  at  a  sample  rate  of  10  Hz.  The  transverse  locations  for  the 
sensors,  measured  outboard  from  the  model  centerline,  were  at  y/B  (distance  outboard 
divided  by  the  transom  beam)  of  0.86,  1.13,  1.37,  1.63,  1.88,  2.12,  and  2.37.  The  average 
height  of  the  sensors  off  the  water  level  was  about  43.5  in  (1 10  em). 
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Figure  5.  Wave  Boom  which  holds  the  ultrasonic  sensors. 


Acoustic  Wave  and  Current  Profiler  (AW AC) 


The  Nortek  Acoustic  Wave  and  Current  Profiler  AW  AC  is  an  acoustic  Doppler 
current  profiler  with  some  added  features.  In  addition  to  the  three  acoustic  beams  angled 
at  25  degrees  from  vertical,  which  are  typically  found  on  an  ADCP,  the  AWAC  system 
has  a  dedicated  vertical  center  beam  which  is  used  to  measure  the  water  surface  through 
Acoustic  Surface  Tracking  (AST).  This  center  beam  transmits  a  short  acoustic  pulse  that 
can  be  finely  resolved,  allowing  for  free  surface  waves  of  short  periods  to  be  accurately 
measured.  The  acoustic  return  at  a  fine  vertical  resolution  may  be  correlated  to  the 
entrained  air  in  the  water.  The  AWAC  is  capable  of  sampling  at  4  Hz  to  capture  the 
surfaee  level;  if  all  bins  are  recorded  to  acquire  acoustic  return  through  the  water  column, 
the  sampling  frequency  is  limited  to  2  Hz.  The  AWAC  (Figure  6)  was  stationary  during 
the  testing,  on  a  bottom  mount  about  halfway  along  the  length  of  the  tank,  located  near 
the  wave  boom  location  whieh  held  the  ultrasonie  sensors.  Measurements  were  made 
over  all  speeds  tested  while  the  AWAC  was  bottom  mounted  under  the  centerline  of  the 
model,  as  well  as  21.75  and  51.75  inches  port  of  eenterline. 
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Figure  6.  Nortek  Acoustic  Wave  and  Current  Profiler  (AW AC)  on  bottom  mount. 


Defoeused  Digital  Imagery  Particle  Image  Veloeimetrv  (DDPIV) 

To  characterize  bubbles,  a  Defoeused  Digital  Particle  Image  Velocimctry 
(DDPIV)  camera  was  employed  in  an  underwater  enclosure.  DDPIV  is  a  volumetric  3D 
measurement  technique  that  is  capable  of  measuring  large  numbers  of  bubbles  at  a  time. 
Illumination  of  the  bubbles  was  provided  by  a  LABest  laser  that  was  synchronized  with 
the  eamera.  Images  from  the  camera  were  recorded  to  hard  drive  on  computers  located  on 
the  eenter  aisle. 

The  DDPIV  eamera  was  placed  near  the  centerline  of  the  basin  on  a  lift  jack 
(Figure  7).  To  measure  bubbles,  the  eamera  needs  the  bubbles  to  be  illuminated  in 
forward  scatter.  To  facilitate  this,  the  eamera  was  placed  at  a  45°  angle,  faeing  upward 
and  towards  the  eenter  aisle  (Figure  8).  The  eamera  was  located  24  in  (0.61  m)  from  the 
centerline  of  the  model,  and  21  in  (0.53  m)  below  the  free  surfaee,  a  depth  that  was 
primarily  designed  to  minimize  the  risk  of  interference  between  the  model  and  the 
eamera.  The  measurement  volume  was  approximately  0.009  in3  (150  mm3),  and  images 
were  collected  at  7  double  frames  per  second.  The  laser  was  mounted  on  the  eenter  aisle, 
with  optics  designed  to  channel  the  light  underwater.  The  beam  was  projected  parallel  to 
the  free  surfaee  at  the  designed  measurement  depth. 
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Figure  7.  DDPIV  camera  on  lift  jack  in  next  to  Carriage  2. 


21"  from  rest  free  surface 

Figure  8.  Plan  view  schematic  of  camera  position  relative  to  transom  model. 


Void  Fraction  Probes 


A  set  of  six  impedance  void  fraction  probes  were  used  to  measure  the  percentage 
of  air  in  the  rooster  tail  generated  by  a  submerged  transom.  The  design  is  based  on  a 
probe  developed  by  Waniewski  (6)  to  measure  void  fraction  in  high-speed,  unsteady, 
multiphase  flows.  Eight  probes  were  built  and  used  by  Coakley  et  al  (7)  in  the  Circulating 
Water  Channel  (CWC)  at  NSWCCD.  In  this  experiment,  only  six  of  the  Coakley  probes 
were  used. 

The  probes  consist  of  two  concentric  electrodes  separated  by  insulation  as  shown 
in  Figure  9.  The  probe  tips  are  aligned  with  the  tow  direction  of  the  model  and  their  small 
dimensions  allow  them  to  respond  to  individual  bubbles.  The  outer  electrode  is 
approximately  0.125  in  (0.3  cm)  in  diameter  and  is  grounded.  A  sinusoidal  voltage  signal 
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of  ±2.5  V  and  an  excitation  of  500  kHz  are  applied  to  the  inner  electrode.  The  impedance 
across  the  two  electrodes  increases  with  void  fraction  (%  of  air)  and  is  mainly  resistive 
for  excitation  frequencies  below  the  megahertz  level.  When  a  bubble  comes  into  contact 
with  the  probe,  the  current  between  the  two  electrodes  decreases  and  voltage  output  of 
the  probe  is  a  large  negative  spike.  The  sampling  rate  of  the  probes  was  set  at  1 8  kHz  and 
was  determined  based  on  limitations  of  the  data  acquisition  system. 


Figure  9.  Cross  seclion  of  1VFM  probe. 


The  probes  were  mounted  on  a  brass  strut  and  aligned  vertically  as  shown  in 
Figure  10.  The  probes  were  positioned  on  the  strut  3.5  in  (9  cm)  apart,  and  the  strut  was 
attached  to  an  aluminum  plate  that  allowed  for  vertical  movement  of  the  strut.  Testing 
conditions  involved  two  vertical  strut  positions.  Data  was  collected  in  vertical  position  1 
and  then  the  strut  was  moved  1.75  in  (4.5  cm)  downward  to  vertical  position  2.  This 
allowed  for  data  to  be  collected  every  1.75  in  (4.5  cm)  for  a  vertical  span  of  19.25  in  (49 
cm).  The  strut  was  attached  to  a  traverse  that  allowed  for  the  longitudinal  position  of  the 
probes  relative  to  the  transom  to  be  controlled.  The  traverse  provided  3  ft  (0.91  m)  of 
longitudinal  movement.  In  addition,  the  traverse  was  mounted  on  slides  attached  to  pieces 
of  aluminum  rail  which  were  part  of  the  structure  used  to  mount  the  void  fraction 
instrumentation  to  the  carriage.  The  slides  were  able  to  move  along  the  rail  providing  a 
full  testing  range  of  8  ft  (2.4  m),  allowing  the  tip  of  the  probes  to  be  positioned  as  close 
as  6  in  (15  em)  and  as  far  as  8.5  ft  (2.6  m)  aft  of  the  transom.  The  longitudinal  locations 
were  chosen  so  that  data  was  collected  throughout  the  white  water  region  of  the  wake  and 
to  ensure  data  was  collected  throughout  the  rooster  tail  from  its  inception  to  its  peak. 
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Figure  10.  Void  fraction  probes  on  brass  strut. 


Figure  1 1  shows  a  schematic  of  the  test  set-up.  The  strut  was  mounted  to  the 
carriage  so  that  the  probes  were  aligned  with  the  tow  direction  of  the  model  and  in  line 
with  the  model  eenterline.  The  probes  were  held  in  one  position  for  each  run,  and  moved 
between  runs.  Figure  1 1  shows  the  ealm  water  waterline  relative  to  vertical  position  1. 
At  this  position  probes  4,  5,  and  7  were  above  the  waterline  at  ealm  water  conditions. 
There  is  no  probe  labeled  number  6  because  the  probe  labeling  system  corresponds  to  the 
eleetronies  card  associated  with  the  probe  and  eleetronies  eard  6  was  not  functioning.  At 
vertical  position  2  the  probe  locations  are  1.75  in  (9  em)  lower  than  vertical  position  1, 
therefore  only  probes  5  and  7  were  above  the  ealm  water  waterline.  Table  4  provides  a 
summary  of  probe  heights  in  inches  relative  to  ealm  water  for  each  vertical  position.  No 
data  was  collected  at  the  probe  height  of  5.75  in  (15  em)  above  the  ealm  water  level 
beeause  Probe  5  stopped  working  while  testing  before  data  at  vertical  position  2  was 
attained.  Void  fraction  data  was  collected  for  one  dry  transom  (7  knots)  and  one  wet 
transom  (8  knots)  condition,  as  shown  in  Table  5. 


Figure  1 1 .  Side  view  schematic  of  test  set-up  (not  to  scale).  The  waterline  shown  represents  probe  positions 
relative  to  calm  water  conditions  at  vertical  position  1. 
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Table  4.  Probe  heights  relative  to  calm  water. 


Height, 

inches 

Probe  # 

Vert  Pos  1 

Vert  Pos  2 

7 

7.5 

5.75* 

5 

4 

2.25 

4 

0.5 

-1.25 

3 

-3 

-4.75 

2 

-6.5 

-8.25 

1 

-10 

-11.25 

*no  data  was  collected  at  this  height 


During  calibrations,  the  gain  for  each  probe  was  set  by  adjusting  the 
potentiometer  in  the  circuit.  The  gain  is  defined  as  the  maximum  value  the  signal  reaehes 
when  the  probe  changes  mediums.  The  gain  for  each  probe  was  set  at  approximately  2V, 
exeept  for  Probe  2.  Probe  2  was  extremely  sensitive  to  gain  ehanges  meaning  that  a 
small  adjustment  to  the  potentiometer  resulted  in  a  large  ehange  in  gain.  Therefore  the 
smallest  gain  achievable  by  Probe  2  was  approximately  4V.  During  testing  it  was 
discovered  that  while  the  gain  had  not  been  physieally  changed  by  adjusting  the 
potentiometer,  the  gains  for  each  probe  were  different  than  the  initial  gains  set  during 
calibrations.  It  is  unclear  why  the  gain  changed,  but  it  is  suspected  that  it  was  a  result  of  a 
flaw  in  the  circuit.  No  attempt  was  made  to  adjust  the  gains  back  to  calibration,  and  it 
was  instead  decided  to  perform  an  extensive  post  calibration  investigation. 


Table  5.  Void  Fraction  Test  Matrix 


Speed 

knts 

Vertical 

Position 

Longitudinal  Position 
inches  aft  of  Transom 

#  of  Runs 

7 

1 

21 

2 

7 

1 

26 

2 

7 

1 

31 

2 

7 

1 

36 

2 

7 

2 

21 

2 

7 

2 

26 

2 

7 

2 

31 

2 

7 

2 

36 

2 

8 

1 

41 

2 

8 

1 

46 

2  !' 

8 

1 

51 

2 

8 

1 

56 

2 

8 

2 

41 

2 

8 

2 

46 

2 

8 

2 

51 

2 

8 

2 

56 

2 

Total  runs  = 

32 
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RESULTS 


Standard  Video  and  Still  Imaging 

Figure  12  shows  still  images  recorded  of  the  transom  wake  during  testing  over  at 
5  knots,  7  knots,  8  knots  and  9  knots.  From  these  images,  it  is  shown  that  the  5  and  7 
knot  conditions  are  wet  transom  conditions,  while  the  8  and  9  knot  conditions  are  dry 
transom  conditions.  There  was  no  significant  change  to  the  free  surface  level  aft  of  the 
transom  at  5  knots,  though  the  disturbance  corresponding  to  the  depth  Froude  number 
change  and  white  water  can  be  seen.  The  development  of  the  “rooster  tail”  aft  of  the  stem 
is  apparent  at  9  knots. 


Figure  12.  Still  images  of  transom  wake  during  testing  over  all  tested  speeds. 


Forces 


For  the  following  figures,  the  forces  and  displacements  were  averaged  over  the 
various  longitudinal  positions  tested  (see  Table  2  above).  Figure  13  shows  the  average 
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drag  force  versus  speed  for  all  of  the  positions  over  the  range  of  speeds  tested.  As 
expeeted,  the  magnitude  of  the  drag  increases  with  increased  speed.  Figure  14  shows  the 
average  vertical  foree  versus  speed  for  all  of  the  positions  over  the  range  of  speeds  tested. 
Vertical  foree  also  increases  with  increased  speed.  Side  foree  was  also  monitored  during 
testing  to  ensure  that  the  model  remained  aligned  properly.  Side  forees  stayed  within  the 
allowable  range  of  1-2%  of  drag  force. 


Model  5673  Drag  Force 


Figure  13.  Drag  force  versus  model  speed  for  all  positions. 


15 


Model  5673  Vertical  Force 


Model  Speed  (knots) 

Figure  14.  Vertical  force  versus  speed  for  all  positions. 


Trim 


Figure  15  shows  the  vertical  displacement  at  the  forward  and  aft  string 
potentiometer  versus  speed  for  all  positions  over  the  range  of  speeds  tested,  where  a 
negative  number  indicates  that  the  bow  has  moved  up.  With  the  bow  moving  up  and  the 
stem  moving  down  as  speed  increase,  the  trim  angle  is  increasing  with  increased  speed. 
Table  6  shows  the  calculated  trim  angle  and  draft  values  at  the  forward  and  aft 
perpendicular. 
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Model  5673  Trim 


Figure  15.  Average  forward  and  aft  trim  values  over  all  speeds. 


Table  6.  Trim  angle  and  draft  at  forward  and  aft  perpendiculars. 


Speed 

ReL 

Fnt 

Trim  Angle 

Tfp  „ 

Tap 

Fn0 

(kts) 

(deg) 

(ft) 

(ft) 

5 

2.41E+07 

0.27 

0.19 

0.97 

1.08 

1.4 

7 

3.37E+07 

0.38 

0.51 

0.93 

1.20 

1.9 

8 

3.85E+07 

0.43 

0.73 

0.90 

1.29 

2.1 

9 

4.33E+07 

0.49 

0.78 

0.90 

1.31 

2.3 

LiDAR 

One  transverse  wave  profile  was  collected  by  the  LiDAR  system  during  each  run. 
The  LiDAR  data  for  the  run,  was  then  averaged  and  eorreeted  for  the  constant  tilt  of  the 
system.  Mean  wake  profiles  of  Model  5673  were  then  generated.  Figure  16,  Figure  17, 
Figure  1 8,  and  Figure  19  show  the  mean  wave  elevations  aft  of  the  model  for  eaeh  of  the 
four  speeds  5,  7,  8  and  9  knots,  respectively.  These  plots  eonfirm  the  expected  increase 
in  wave  elevation  with  increase  in  speed.  The  data  from  these  runs  were  also  used  to 
create  two-dimensional  and  contour  plots  in  an  effort  to  visually  re-create  the  wake  of 
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Model  5673  and  arc  shown  in  Figure  20,  Figure  21,  Figure  22,  and  Figure  23  for  each  of 
the  four  speeds  5,  7,  8  and  9  knots,  respectively.  These  plots  look  similar  to  the  still 
images  shown  in  Figure  12,  with  very  little  change  in  the  free  surface  at  the  5  knot 
condition,  and  the  development  of  the  “rooster  tail”  at  the  9  knot  condition. 


Wake  Profiles;  5  kts  Wet  Transom,  Single  Runs 


Figure  16.  Mean  wake  Profile  of  model  5673,  5  knots. 


Wake  Profiles.  7  kts  Wet  Transom,  Single  Runs 


Figure  17.  Mean  wake  profile  of  model  5673,  7  knots. 
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Height  (m)  Height  (m) 


Wake  Profiles  8  kts  Dry  Transom,  Single  Runs 


Wake  Profiles  9  kts,  Dry  Transom,  Single  Runs 


Figure  19.  Mean  wake  profile  of  model  5673,  9  knots. 
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Figure  20.  2-D  characterization  of  model  5673  wake,  5  knots. 


Contour  Plot,  7  knots,  Wet  Transom 


sttxl  or  lidar  (m) 

Figure  21.  2-D  characterization  of  model  5673  wake,  7  knots. 
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Contour  Plot,  8  knots,  Dry  Transom 


stbd  of  Ildar  (m) 


Figure  22:  2-D  characterization  of  model  5673  wake,  8  knots 


Contour  Plot,  9  knots,  Dry  Transom 


Figure  23:  2-D  characterization  of  model  5673  wake,  9  knots. 


DLP-Enhaneed  QViz 


Analysis  of  the  results  from  the  DLP-enhaneed  QViz  system  is  ongoing. 
Fluctuations  in  the  free  surface  profiles,  as  well  as  averaged  free  surface  contours,  arc 
currently  being  studied.  The  ongoing  data  analysis  involves  edge  detection  algorithm 
development,  to  identify  the  free  surface  in  each  of  the  raw  images,  as  well  as  calibration 
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method  advancement.  These  results  will  allow  comparisons  to  be  made  with  other 
measurements  made  during  testing,  as  well  as  with  full-scale  field  measurements  of  the 
transom  wave  of  the  R/V  Athena  reported  in  Fu,  et  al  (8  and  9). 

Wavecuts  (Senix  Ultrasonic  Sensors) 

The  data  from  the  Senix  ultrasonic  sensors  was  filtered  to  remove  dropouts  and 
obtain  a  smooth  wave  record.  Figure  24  shows  an  example  of  the  data  that  was  collected 
for  the  9  knot  condition.  The  black  line  represents  the  original  signal,  while  the  red  line 
shows  the  smoothed,  filtered  data.  This  data  was  then  used  to  calculate  the  wavemaking 
resistance  coefficient,  results  of  which  arc  tabulated  in  Table  7  as  Cw,  along  with  the 
other  resistance  coefficients.  These  coefficients  are  plotted  in  Figure  25.  The  wave 
resistance  coefficient  tends  to  increase  as  speed  increases  for  this  transom  model.  Also 
included  in  this  table  are  the  total  resistance  coefficients  (Ct),  the  frictional  resistance 
coefficients  (Cf),  and  the  residuary  resistance  coefficients  (Cr).  The  total  resistance 
coefficient  is  calculated  using  the  measured  drag  force  and  the  wetted  surface  area 
(calculated  from  the  measured  sinkage  and  trim).  The  frictional  resistance  coefficient  is 
calculated  using  the  ITTC  57  formula  (PNA,  10).  Residuary  resistance  is  calculated  by 
subtracting  the  frictional  resistance  coefficient  from  the  total  resistance  coefficient.  The 
wavemaking  resistance  coefficient  makes  up  part  of  the  residual  resistance,  along  with 
the  eddy  resistance.  All  resistance  coefficients  are  non-dimensionalized  by  the  static 
wetted  surface  area  at  the  1  foot  draft,  zero  trim  condition  (132.36  ft2,  12.3  m2). 


Figure  24.  Example  of  original  and  ultrasonic  wave  record  for  9  knot  condition.  Sonic  #1  is  closest  to 

centerline  of  model. 
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Table  7.  Table  of  resistance  coefficients. 


Speed 

Fnt 

R« 

Ct 

Cf 

Cr 

Cw 

(kts) 

non-dim 

(lbs) 

5 

0.27 

95.9 

0.0105 

0.0026 

0.0079 

0.00001 

7 

0.38 

217.3 

0.0121 

0.0025 

0.0097 

0.00203 

8 

0.43 

266.9 

0.0114 

0.0024 

0.0090 

0.00278 

9 

0.49 

295.5 

0.0100 

0.0024 

0.0076 

0.00340 

Figure  25.  Coefficients  of  resistance  for  transom  model  computed  across  range  of  FToude  numbers  tested. 


Acoustic  Wave  and  Current  Profiler  (AWAC) 


Figure  26  shows  the  water  level  data  collected  by  the  AWAC  at  the  model 
eenterline  position,  with  time  converted  to  distance  from  the  bow.  Each  panel  shows  2 
separate  runs  for  eaeh  speed,  with  the  5  knot  ease  in  the  top  panel  and  the  9  knot  ease  in 
the  bottom  panel.  The  aft  edge  of  the  model  can  be  seen  in  the  plot  at  approximately  30 
feet  (9.1  m).  The  wake  behind  the  ship  is  difficult  to  resolve  from  the  plots,  likely  due  to 
the  sampling  rate  of  4  Hz.  The  transverse  wave  is  apparent  in  the  7,  8,  and  9  knot  plots, 
and  the  length  of  this  wave  increases  with  vessel  speed,  as  expeeted. 
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Figure  26.  Water  level  time  series 


Figure  27  shows  the  aeoustie  return  for  the  ecnter  beam  of  the  AWAC  in  counts 
while  the  instrument  was  in  the  eenterline  position.  Eaeh  figure  shows  a  single  run,  with 
the  time  converted  to  distance  past  the  bow.  As  with  Figure  26,  the  aft  edge  of  the  model 
can  be  seen  in  the  plot  at  approximately  30  feet  (9.1  m).  Eaeh  panel  shows  a  different 
speed,  with  5  knots  in  the  top  panel  and  9  knots  in  the  bottom  panel.  The  blaek  line 
shows  the  level  of  maximum  acoustic  return,  which  is  the  water  surface  or  the  outline  of 
the  model.  From  these  plots,  it  can  be  seen  that  the  model  sinks  lower  at  greater  speeds, 
which  was  verified  through  the  draft  measurements  taken  during  testing.  Assuming 
aeoustie  return  is  related  to  bubble  density,  it  appears  that  more  bubbles  are  present  lower 
in  the  water  eolumn  at  lower  speeds  than  at  greater  speeds. 

These  plots  can  be  compared  with  the  still  images  of  the  flow  aft  of  the  transom 
(Figure  12),  which  show  that  at  5  knots,  the  transom  flow  is  relatively  flat  with  many 
bubbles  behind  the  hull.  From  the  still  photos,  it  appears  that  as  the  speed  increases,  the 
density  of  bubbles  observed  decreases,  resulting  in  a  eleaner  flow  with  bubbles  only  at 
the  surface,  which  tends  to  agree  with  what  is  seen  in  the  acoustic  return  plots.  The 
aeoustie  plots  eontain  data  in  raw  counts,  so  further  analysis  is  necessary  to  facilitate  the 
direct  comparison  to  bubble  size  and  density. 
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distance  (ft) 

Figure  27.  AWAC  return  from  centerline  for  5,  7,  8,  and  9  knots. 


Figure  28  and  Figure  29  show  the  acoustic  return  for  two  other  locations,  21.75  in 
(55.2  cm)  port  of  the  model  centerline  and  51.75  in  (131.4  cm)  port  of  the  model 
centerline.  The  AWAC  return  at  the  21.75”  port  of  centerline  position  looks  similar  to 
the  return  at  the  centerline  position,  with  a  higher  intensity  return  penetrating  deeper  at 
the  lower  speeds  than  at  the  greater  speeds.  The  overall  acoustic  return  is  much  lower  at 
the  51.75  in  (131.4  cm)  port  of  centerline  position  than  at  the  centerline  or  21.75  in  (55.2 
cm)  port  of  centerline  position,  with  very  little  return  other  than  the  water  surface, 
particularly  for  the  higher  speeds.  This  trend  makes  sense  when  examining  the  still 
photos  of  the  flow  (Figure  12),  where  the  bubbles  on  the  surface  spread  less  across  the 
beam  of  the  wake  with  greater  speeds. 


25 


9kts- depth  (ft)  8kts- depth  (It)  7KtS- depth  (tt)  Skts- depth  (ft)  9Ms- depth  (ft)  flkts- depth  (ft)  7Ws- depth  (ft)  5kts- depth  (ft) 


distance  (ft) 

Figure  28.  AW  AC  return  from  21.75  in  (55.2  cm)  port  of  centerline  for  5,7,8,  and  9  knots. 


distance  (ft) 

Figure  29.  AWAC  return  from  51.75  in  (131.4  cm)  port  of  centerline  for  5,7,8,  and  9  knots. 
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Defocuscd  Digital  Imagery  Particle  Image  Velocimctrv  (DDPIV) 


Because  the  camera  was  fixed  to  the  basin,  the  DDPIV  test  procedure  was  to  take 
images  continuously  during  the  carriage  transit,  and  then  select  the  frames  which  had 
bubbles  present  for  processing.  In  most  cases,  the  camera  would  only  capture  one  frame 
that  had  bubbles  present  for  each  run.  Although  the  carriage  was  run  at  speeds  from  7  to 
9  knots,  bubbles  were  only  present  in  the  imaging  volume  at  the  higher  speeds  (8  and  9 
knots).  It  should  be  noted  that  the  camera  was  located  24  in  (0.61  m)  from  the  centerline 
of  the  model  and  21  in  (0.53  m)  below  the  surface  so  as  to  avoid  interference  with  the 
model,  as  indicated  in  the  previous  section.  The  contour  plots  of  the  raw  signal  from  the 
AW  AC  shown  above  indicate  that  there  are  more  bubbles  present  in  this  region  at  lower 
speeds  than  higher  speeds,  which  contradicts  these  results.  It  may  be  that  though  there 
are  fewer  bubbles,  they  arc  also  smaller  and  more  persistent  at  the  higher  speeds, 
allowing  the  camera  to  capture  these  images.  The  larger  bubbles  produced  at  lower 
speeds  will  surface  more  quickly.  Also,  the  AWAC  data  is  in  raw  form,  and  needs 
further  processing  before  it  can  be  used  as  a  proxy  for  void  fraction.  Nine  cases  were 
recorded  at  8  knots  and  eight  cases  at  9  knots.  A  sample  of  the  raw  bubble  image  is 
shown  in  Figure  30. 


.  «.  *  •  *•  «  4b 

■  - .  .  •  ,  #  ♦  ‘  *  ****.. 

■  *  \ '  *  *  *-»  •  • 

•  i. 

#  k  ~  .  .• 

•Vts/  -..V 

TMKt : 


Figure  30.  Sample  of  raw  bubble  image  from  DDPIV. 

The  bubble  size  distribution  is  shown  in  Figure  3 1 .  In  both  the  8  knot  and  9  knot 
cases,  the  peak  of  the  distribution  is  around  200  microns.  The  rise  in  the  distribution  at 
the  small  diameters  is  most  likely  due  to  spurious  data.  Although  the  distribution  appears 
to  be  slightly  shifted  towards  smaller  size  at  the  higher  speed,  it  is  unlikely  to  be  a 
significant  difference.  The  bubble  populations  were  then  averaged  in  columns  parallel  to 
the  free  surface  and  integrated  to  obtain  void  fraction.  The  void  fraction  distribution  for 
the  9  knot  case  is  shown  in  Figure  32.  The  concentration  of  high  void  fraction  at  depth  is 
most  likely  indicative  of  air  entrained  at  the  toe  of  the  stem  wave  which  was  close  to  the 
measurement  volume.  Corresponding  velocity  data  could  confirm  this  hypothesis. 
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Figure  3 1 .  Bubble  size  distribution  from  DDPIV  for  8  and  9  knot  cases. 
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Figure  32.  Void  fraction  distribution  from  DDPIV  for  knot  case. 


Void  Fraction  Probes 


During  the  experiment  probes  began  a  run  either  in  air  or  water,  and  depending  on 
the  test  condition  some  probes  remained  out  of  the  water  for  the  entire  run.  All  probes 
were  zeroed  in  water  regardless  of  whether  the  probe  was  submerged  or  not.  For  those 
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probes  above  the  calm  water  waterline,  a  cup  of  water  was  used  to  submerge  the  probe 
while  zeroing,  since  the  baseline  condition  was  considered  as  submerged. 

A  voltage  signal  was  observed  in  the  data  that  would  switch  between  two 
extremes,  which  represents  a  phase  change.  A  probe  that  was  underwater  during  a  run 
would  primarily  be  at  the  upper  extreme  and  would  spike  to  the  lower  extreme  when  the 
probe  encountered  a  bubble,  as  shown  in  Figure  33  and  Figure  34.  A  probe  that  was  out 
of  the  water  during  a  run  would  primarily  be  at  the  lower  extreme  and  would  spike  to  the 
upper  extreme  when  the  probe  encounter  spray  as  show  in  Figure  35.  Figure  36  shows 
data  collected  for  a  probe  that  was  near  the  free  surface. 

To  determine  the  void  fraction  level  from  the  data  a  voltage  threshold  was 
applied.  All  points  with  a  value  smaller  than  the  threshold  (i.c.  lie  below  the  threshold) 
represent  time  that  the  probe  is  in  air.  Data  points  with  a  value  larger  than  the  threshold 
(i.e.  lie  above  the  threshold)  represent  time  the  probe  is  in  water.  Therefore,  the  void 
fraction  level,  represented  by  a  percentage,  was  determined  by  the  ratio  of  the  number  of 
data  points  that  are  below  the  threshold  (time  the  probe  is  in  air)  to  the  total  number  of 
data  points  collected.  Each  probe  had  its  own  unique  threshold  that  was  determined 
during  calibration.  Details  of  the  probe  calibration  can  be  found  in  Appendix  A. 

Before  void  fraction  levels  were  calculated  the  median  value  of  each  signal  was 
removed  to  ensure  that  all  values  were  properly  zeroed  before  analysis.  The  probes  and 
electronics  were  sensitive  to  changes  in  the  test  environment  and  although  zeros  were 
collected  before  each  run  occasionally  the  slightest  change  in  the  system  would  cause  the 
signal  to  shift.  Removing  the  median  ensured  the  analyzed  signal  was  always  zeroed  at 
zero  before  analysis.  In  order  for  the  threshold  method  to  be  applied  to  the  signal  it  was 
imperative  that  the  signal  was  properly  zeroed  so  that  the  signal  noise  fluctuated  about 
zero. 


Probe  2:  Bknte,  x=56",  vertical  pos  t 


Time,  seconds 


Figure  33.  Example  of  data  for  an  underwater  probe  that  sees  air  bubbles.  The  red  line 
represents  the  threshold  used  in  analysis. 
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Probe  2:  8knts,  x=56\  vertical  pos  1 


Figure  34.  A  close  up  of  data  presented  in  Figure  33.  The  circles  represent  the  data  points 
collected.  The  red  line  represents  the  threshold  used  in  analysis. 


Figure  35.  Example  of  data  for  a  probe  that  is  in  air  and  sees  water  droplets.  The  red  line 

represents  the  threshold  used  in  analysis. 
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Probe  5:  Bknts,  x=51 ",  vertical  pos  1 
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Time,  seconds 


Figure  36.  Example  of  data  for  a  probe  that  moves  in  and  out  of  the  free  surface.  The  red 
line  represents  the  threshold  used  in  analysis. 

Void  fraction  data  was  processed  for  all  runs,  and  then  averaged  over  the 
duplicate  runs,  which  is  the  data  presented  in  this  section.  Figure  37  and  Figure  38  show 
the  variation  of  void  fraction  with  probe  height  and  longitudinal  distance  aft  of  the 
transom  for  the  two  tested  speeds  of  7  and  8  knots.  The  probe  heights  shown  in  the  y-axis 
arc  relative  to  calm  water  and  the  legend  refers  to  the  distance  in  inches  aft  of  the 
transom,  with  a  larger  number  indicating  a  greater  distance  aft  of  the  transom.  Error  bars 
are  provided  for  each  data  point  and  were  calculated  using  a  95%  confidence  interval  and 
assuming  a  t-distribution.  Details  of  this  error  analysis  can  be  found  in  Appendix  A.  The 
majority  of  the  uncertainty  is  found  in  the  mid  to  upper  range  of  void  fraction  levels, 
specifically  50-90%.  This  is  expected  as  both  Figure  37  and  Figure  38  show  that  void 
fraction  increases  rapidly  through  this  region  showing  that  a  small  change  in  probe  height 
results  in  a  large  change  in  percentage  of  void  fraction. 
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Probe  Height  Relative  to  Calm  Water,  inches  Prob®  Relative  to  Calm  Water,  inches 


7  knts:  Wet  Transom 
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Figure  37.  Probe  Height  vs  Void  Fraction  for  7  knots. 


8  knts:  Dry  Transom 
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Figure  38.  Probe  Height  vs.  Void  Fraction  for  8  knots. 
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Figure  39  and  Figure  40  show  the  void  fraction  eontour  plots  for  both  the  wet  and 
dry  transom  conditions.  The  contour  lines  represent  the  pereent  of  air  in  the  flow;  the 
50%  contour  line  is  a  rough  approximation  of  the  free  surfaee.  The  spaeing  of  the  eontour 
lines  was  chosen  based  on  the  error  analysis.  The  spacing  is  equal  to  the  largest  range  of 
error  for  a  given  region.  For  example  for  void  fraction  levels  ranging  from  1-10%  the 
largest  error  calculated  in  this  range  was  2%,  therefore,  between  1-10%  void  fraction  the 
contour  lines  shown  arc  for  every  2%.  For  void  fraction  levels  ranging  from  10-30%  the 
spacing  is  5%  and  from  30-100%  the  spaeing  is  20%.  The  x-axis  represents  the 
longitudinal  distance  relative  to  the  transom  with  negative  values  indicating  a  distance  aft 
of  the  transom.  The  y-axis  provides  the  probe  height  relative  to  calm  water.  All  data  was 
taken  along  the  model  centerline  and  there  was  no  variation  in  transverse  location.  The 
black  dots  represent  the  probe  locations  where  data  was  collected.  Longitudinal  data  was 
taken  every  5  inches  and  vertical  data  every  1.75  inches. 

The  contour  plots  concur  with  the  visual  observations  made  from  the  photographs 
in  Figure  12.  For  the  7  knot,  wet  transom  case,  the  linear  trend  of  the  free  surfaee  is 
confirmed.  The  plot  in  Figure  39  shows  that  the  free  surface  gradually  spans  a  height  of 
approximately  5  inches  over  a  15  inch  longitudinal  distance.  In  contrast,  the  dry  transom 
case  at  8  knots  has  a  steeper  climb  spanning  5.25  inches  vertically  over  a  9  ineh 
longitudinal  distance.  The  plot  in  Figure  39  indicates  the  presence  of  a  small  vortex  near 
x=-31”  and  y=-10”  which  corresponds  to  the  knob  seen  in  Figure  37  at  x=-31”  and  y=- 
10”.  This  vortex  is  not  present  at  the  8  knot,  dry  transom  condition.  While  this  data  shows 
the  presence  of  a  vortex,  it  is  difficult  to  draw  conclusions.  It  is  important  to  remember 
the  resolution  of  the  grid  used  in  analysis  is  rough  and  each  condition  was  only  repeated 
onee.  This  is  an  area  to  be  explored  in  more  detail  in  the  future. 
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Figure  39.  Void  fraction  contours  at  wet  transom  condition,  7  knots. 


8  knts:  Dry  Transom 


Figure  40.  Void  fraction  contours  at  dry  transom  condition,  8  knots. 
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CONCLUSIONS 


A  laboratory  model  capable  of  producing  full-scale  transom  breaking  waves 
similar  in  behavior  to  those  of  full-scale  naval  combatants  has  been  designed,  fabricated 
and  tested.  Multiple  measurement  methods  have  been  implemented  and  presented,  and 
comparisons  between  the  still  images  of  the  wake  and  the  LiDAR  data,  AWAC  data  and 
void  fraction  data  seem  to  agree  well.  The  initial  experimental  work  to  document  and 
characterize  a  large  breaking  transom  wave  in  calm  water  over  a  range  of  transom  depth 
Froudc  numbers  has  been  completed,  and  analysis  of  data  continues  to  provide  further 
comparisons  between  measurements,  as  well  as  insight  into  the  physics  of  the  transom 
stem  wave. 
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APPENDIX  A:  VOID  FRACTION  PROBE  CALIBRATION  PROCEDURE  AND 

ERROR  ANALYSIS 


Set-Up 


The  probes  were  calibrated  using  the  set  up  shown  in  Figure  A-l.  A  vertical 
column  was  filled  with  water  and  had  two  pressure  taps  loeated  12  in  (30.5  cm)  apart  that 
were  eaeh  attached  to  a  differential  pressure  gage.  Air  stones  connected  to  an  air  line  that 
was  connected  to  an  air  compressor  were  used  to  injeet  bubbles  into  the  eolumn  of  water, 
creating  a  two-phase  flow.  The  air  stones  were  plaeed  at  the  bottom  of  a  half  cylinder 
piece  of  clear  PVC  tubing  that  was  mounted  flush  against  one  wall  of  the  tank.  This 
allowed  the  bubbly  flow  to  be  confined  to  the  region  inside  the  tube  as  shown  in  Figure 
A-2.  The  probe  was  inserted  in  the  top  of  the  half  cylinder  and  aligned  with  the  flow.  Due 
to  probe  geometry  the  probe  could  only  be  placed  a  few  inehes  below  the  uppermost 
pressure  tap.  The  air  compressor  was  set  at  a  constant  pressure  and  the  pressure  in  the  test 
cylinder  was  controlled  by  the  valve  shown  in  Figure  A-3.  Every  quarter  turn  of  the  valve 
was  marked  and  data  was  eolleeted  through  one  full  turn  of  the  knob,  thus  eight  valve 
locations. 


pressure  taps 


Side  View 


Front  View 


Figure  A-l.  Water  column  used  to  calibrate  1VFM. 
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Figure  A-2.  Pictures  showing  region  of  bubbly  flow. 


Figure  A-3.  Regulator  valve. 


Procedure 

The  probes  were  calibrated  using  a  threshold  method  used  by  Waniewski  (6). 
Denoting  the  distance  between  the  pressure  taps  as  H  and  the  differential  pressure 
measured  by  the  pressure  gages  as  h,  the  steady  state  void  fraction,  a,  was  determined 
using  Equation  A-l. 


(A-l) 
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A  threshold  value  was  chosen  so  that  the  percentage  of  void  fraction  data  that  was  below 
this  threshold  value  was  equal  to  the  percentage  void  fraction  (a)  ealeulated  with  the 
above  equation. 

Figure  A-4  provides  an  example  of  calibration  data  for  one  of  the  void  fraction 
probes.  The  void  fraction  measured  by  the  probes  using  the  threshold  method  is  shown  in 
the  y  axis  versus  the  void  fraction  ealeulated  using  the  data  collected  from  the  pressure 
gages.  The  threshold  values  were  varied  to  determine  whieh  value  provided  the  elosest  fit 
to  the  ideal  line  represented  in  the  graph  by  the  solid  blue  line.  It  was  determined  for  this 
probe  that  a  threshold  of  -0.8V  was  the  best  mateh.  Eaeh  probe  was  individually 
calibrated  and  its  corresponding  threshold  was  determined. 


Probe  3:  Test  22 


Figure  A-4.  Example  of  void  fraction  probe  calibration  data. 


Error  Analysis 

Before  calibrations  the  gain  for  eaeh  probe  was  set  by  adjusting  the  potentiometer 
in  the  circuit.  To  set  the  gain  eaeh  probe  was  zeroed  while  no  bubbles  were  present  in  the 
flow.  Bubbles  were  then  introduced  into  the  flow  so  that  the  probe  signal  contained  the 
spikes  eharaeteristie  of  the  presence  of  air.  The  gain  was  then  adjusted  so  that  the  spikes 
peaked  at  a  desired  voltage.  The  largest  spike  for  eaeh  probe  was  set  at  approximately 
2V,  except  for  Probe  2.  Probe  2  was  very  sensitive  to  gain  ehanges  meaning  that  a  small 
turn  of  the  potentiometer  resulted  in  a  large  change  in  gain.  Therefore  the  smallest  gain 
achievable  by  Probe  2  was  approximately  4V. 

Onee  testing  began  it  was  discovered  that  while  the  gain  had  not  physieally  been 
ehanged  by  adjusting  the  potentiometer  the  gains  for  eaeh  probe  were  different  than  the 
initial  gains  set  during  calibrations.  Table  A-l  summarizes  the  gains  for  eaeh  probe 
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during  calibrations  and  testing.  It  is  unclear  how  the  gain  changed,  but  it  is  suspected  that 
it  has  to  do  with  an  error  in  the  circuit.  The  circuit  appears  to  be  very  sensitive  to  changes 
in  conductivity  and  resistance  and  by  changing  environments  from  the  calibration  lab  to 
the  carriage  the  circuit  reacted  to  the  new  environment  and  something  caused  the  gains  to 
change.  No  attempt  was  made  to  adjust  the  gains  back  to  calibration  levels  while  testing 
for  concern  that  any  changes  would  add  further  uncertainty.  Instead  a  decision  was  made 
to  perform  an  extensive  post  calibration  investigation. 


Table  A-l .  Summary  of  gains  during  calibrations  and  testing. 


Probe 

Initial  Calibration 

Testing 

Post  Cal  1 

Post  Cal  2 

Gain  (V) 

Gain  (V) 

Gain  (V) 

Gain  (V) 

1 

2.0 

1.5 

1.8 

2.0 

2 

4.3 

6.0 

4.0* 

4.3 

3 

1.95 

1.2 

1.35 

2.0 

4 

2.15 

1.4 

1.65 

2.15 

5 

2.05 

5.0 

** 

7 

2.15 

1.8 

1.9 

2.0 

*  Probe  2  had  a  gain  higher  than  the  limit  (10V),  it  was  manually  set  to  4V 
before  data  for  Post  Calibration  1  was  collected 

**Probe  5  was  damaged  during  testing  and  has  not  been  repaired. _ 


After  testing  was  completed  the  instrumentation  was  returned  to  the  calibration 
lab  and  set  up  to  repeat  the  initial  calibrations  performed  before  the  test.  Post  Calibration 
1  was  performed  “as  is”  with  no  adjustments  made  to  the  gain  allowing  the 
instrumentation  to  be  as  close  to  carriage  conditions  as  possible,  however  Table  A-l 
shows  that  the  gains  were  still  different  than  testing  conditions.  The  gain  for  Probe  2  was 
so  large  that  it  surpassed  the  limit  of  the  electronics  and  the  signal  was  clipped. 
Therefore,  a  post  cal  could  not  be  performed  on  the  probe  as  is.  Instead  the  gain  was  set 
to  4V  for  Post  Calibration  1 .  Post  Calibration  2  was  performed  by  adjusting  the  gain  as 
best  as  possible  back  to  the  values  present  in  the  initial  calibration.  For  each  set  of  gains 
tested  a  corresponding  set  of  threshold  values  were  determined.  These  values  are 
summarized  in  Table  A-2.  The  results  of  the  post  calibrations  show  that  threshold  varies 
with  gain.  However,  the  lack  of  a  clear  trend  between  the  two  indicates  that  threshold  is 
not  solely  dependent  on  gain  and  there  arc  other  factors  that  may  affect  the  threshold 
values  chosen. 
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Table  A-2.  Summary  of  calibrations,  gains,  and  threshold  values. 


Probe 

Calibration 

Gain,  V 

Threshold,  V 

1 

Initial 

2 

0.6 

Post  1 

1.8 

0.55 

Post  2 

2 

0.55 

2 

Initial 

4.3 

1.5 

Post  1 

4 

1.3 

Post  2 

4.3 

1.2 

3 

Initial 

1.95 

0.8 

Post  1 

1.35 

0.4 

Post  2 

2 

0.4 

4 

Initial 

2.15 

0.8 

Post  1 

1.65 

0.45 

Post  2 

2.15 

0.55 

7 

Initial 

2 

0.9 

Post  1 

1 

0.4 

Post  2 

2 

0.5 

The  test  data  was  processed  using  the  three  different  sets  of  threshold  values 
determined  from  the  initial  and  post  calibrations.  The  void  fraction  levels  at  each 
longitudinal  location  and  probe  heights  were  compared  and  are  shown  in  Figures  A-5 
through  A- 12.  For  both  speeds  at  all  longitudinal  locations,  the  trend  of  the  data  is 
consistent  regardless  of  which  set  of  threshold  values  arc  used.  The  variation  in  threshold 
values  simply  shifts  the  data.  There  is  less  variation  among  the  different  threshold  values 
for  the  wet  transom  ease  at  7  knots  than  there  is  at  the  dry  transom  ease  at  8  knots.  Since 
the  trend  of  the  data  was  the  same  regardless  of  threshold,  the  average  void  fraction 
profile  at  each  longitudinal  location  was  calculated  and  that  is  the  data  presented  in  this 
report. 

The  random  uncertainty,  P,  was  calculated  at  each  probe  location  for  all 
longitudinal  locations  using  equation  A-2. 

/>=/<>,  (A-2) 

The  value  for  t  is  determined  assuming  a  95%  confidence  interval  and  N-l  degrees  of 
freedom.  N  represents  the  sample  size  and  in  this  ease  N=3  since  there  are  3  different 
threshold  values  for  each  probe  that  were  used  to  process  the  data.  Sj  is  the  sample 
standard  deviation  and  is  determined  by  equation  A-3. 


5,  = 


N- 1 


i=\ 


(A-3) 


The  random  uncertainty  was  calculated  for  both  speeds  at  each  longitudinal  location  and 
probe  height  and  is  summarized  in  Table  A-3.  These  values  were  used  for  the  error  bars 
shown  in  Figure  37  and  Figure  38. 
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Table  A-3.  Summary  of  Random  Uncertainty  (P)  values. 


Probe 

Probe  Height 
relative  cairn  water,  in 

8  knots 

7  knots 

x=56 

x=51 

x=46 

x=41 

x=36 

x=31 

x=26 

x=21 

1 

-11.75 

-10 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0005 

0.0063 

0.0060 

0.3063 

0.0165 

0.0452 

0.0312 

0.0612 

2 

-8.25 

-6.5 

0.0001 

0.0021 

0.0005 

0.0189 

0.0052 

0.1402 

0.0414 

0.6939 

0.0244 

0.0605 

0.0746 

0.5621 

0.1774 

0.6831 

0.5995 

1.3411 

3 

-4.75 

-3 

0,0753 

0.2193 

0.2983 

0.6255 

1.3189 

1.1848 

3.0530 

4.4714 

0.1319 

0.2649 

0.6034 

4.4489 

4.2405 

20.4126 

15.5476 

13.3245 

4 

-1.25 

0.5 

0.3916 

0.5269 

0.7534 

1.0242 

1.2708 

3.2996 

16.9517 

6.9186 

4.0106 

7.1747 

7.2085 

5.8156 

5.2721 

3.4388 

2.8431 

2.2625 

7 

2.25 

4 

1.4404 

2.5598 

3.4408 

7.6642 

12.5197 

5.8957 

3.2675 

0.1123 

4.4715 

1.9652 

9.4828 

17.5811 

4.7959 

7.9311 

5.8989 

0.8602 
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Calibration  Comparison  Plots 


8  knts  x=56" 


Figure  A-5.  Probe  Height  vs.  Void  Fraction  at  8  knots  for  longitudinal  location  x=56” 


8kntsx=51" 


Figure  A-6.  Probe  Height  vs.  Void  Fraction  at  8  knots  for  longitudinal  location  x=51” 
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8  krtfs  x=46" 


Figure  A-7.  Probe  Height  vs.  Void  Fraction  at  8  knots  for  longitudinal  location  x=46” 


8  knts  x=41" 


Figure  A-8.  Probe  Height  vs.  Void  Fraction  at  8  knots  for  longitudinal  location  x=41” 
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7  knts  x=36" 


Figure  A-9.  Probe  Height  vs.  Void  Fraction  at  7  knots  for  longitudinal  location  x=36”. 


7  knts  x=31" 


Figure  A- 10.  Probe  Height  vs.  Void  Fraction  at  7  knots  for  longitudinal  location  x=31”. 
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7  knts  x=26" 


Figure  A-l  1.  Probe  Height  vs.  Void  Fraction  at  7  knots  for  longitudinal  location  x=26” 


7kntsx=21" 


Figure  A-l 2.  Probe  Height  vs.  Void  Fraction  at  7  knots  for  longitudinal  location  x=21” 
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